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Processes involving alpha particles and alpha-like nuclei comprise a major part of stellar nucleosynthesis and 
hypothesized mechanisms for thermonuclear supemovae. In an effort towards understanding alpha processes 
from first principles, we describe in this letter the first ab initio calculation of alpha-alpha scattering. We use 
lattice effective field theory to describe the low-energy interactions of nucléons and apply a technique called 
the adiabatic projection method to reduce the eight-body System to an effective two-cluster System. We find 
good agreement between lattice results and experimental phase shifts for S-wave and D-wave scattering. The 
computational scaling with particle number suggests that alpha processes involving heavier nuclei are also 
within reach in the near future. 

PACS numbers: 21.10.Dr, 21.30.-x, 21.60.De 
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INTRODUCTION 

In recent years there has been much progress in ab initia nuclear scattering and reactions involving light nuclei iHHa as well 
as medium-mass nuclei naiïi. However, the computational scaling for most numerical methods increases significantly when the 
projectile nucléus has more than a few nucléons. Therefore it remains a challenge to address important alpha processes relevant 
for stellar astrophysics such as ^He -f '^He scattering and -I- "^He scattering and radiative capture, as well as carbon and 
oxygen burning in massive star évolution and in thermonuclear supernovae ÜHÏôl. 

In this Work we describe lattice calculations for which the computational scaling of the Ai-body + A 2 -boAy problem is 
roughly {Ai -b ^ 2 )^, mild enough to make first principles calculations of alpha processes possible. We use the formalism of 
lattice effective field theory miHii and make use of a technique for elastic scattering and inelastic reactions on the lattice called 
the adiabatic projection method OMIll. In the following we présent the first ab initia calculation of '^He -b '^He scattering, 
going up to next-to-next-to-leading order in chiral effective field theory. See Ref. ll^ for a recent review of chiral effective field 
theory. We find good agreement with experimental data for the S'-wave and H-wave phase shifts GIMl. 


ADIABATIC PROJECTION METHOD 

The adiabatic projection method addresses the cluster-cluster scattering problem on the lattice by using Euclidean time pro¬ 
jection to construct an effective two-cluster Hamiltonian for the participating clusters. By Euclidean time projection, we mean 
multiplication by exp(—fïr), where H is the underlying microscopie Hamiltonian. For simplicity of explanation we refer to a 
continuons Euclidean time parameter r, even though the actual lattice calculations use discrète time steps. 

Our starting point is an periodic spatial lattice. We take a set of initial two-alpha States \R) labeled by their séparation 
vector R, as illustrated in Fig.lu We take the initial alpha wave functions to be Gaussian wave packets so that at large séparations 
they factorize as a tensor product of two individual alpha clusters, 

( 1 ) 

r 

The summation over r produces two-alpha States with total momentum equal to zéro. Rather than dealing with a large array of 
three-dimensional vectors R, we find it convenient to project onto spherical harmonies with angular momentum quantum 

numbers L, Lz, 


\R)YI^^ ^^YL,LAR')Sn.\R'\\R')- (2) 

R' 

We then use Euclidean time projection to form dressed cluster States, 

=exp(-iÏT)|iî)^’^^ (3) 


We are using units where H and c are set to 1. The évolution in Euclidean time automatically incorporâtes the induced deformation 
and polarization of the alpha clusters as they corne near each other. 

With these dressed cluster States, we compute matrix éléments of the full microscopie Hamiltonian with respect to the dressed 
cluster States, 


Since the dressed cluster States are not orthogonal, we construct a norm matrix, 

The radial adiabatic Hamiltonian is defined as the matrix product. 


(4) 


(5) 


[h: 


t\r,R' 


Nr ^HrN~ 


1 L,Lz 


J R,R' 


(6) 


In the limit of large projection time r, the spectrum of the adiabatic Hamiltonian reproduces the low-energy finite volume 
spectrum of the microscropic Hamiltonian H. In Ref. El it was shown that in the asymptotic région where the alpha clusters 
are widely separated, the adiabatic Hamiltonian reduces to a simple two-cluster Hamiltonian with only infinite-range interactions 
such as the Coulomb interaction between the otherwise non-interacting clusters. 
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FIG. 1. Initial State |iî) of two alpha particles separated by the displacement vector R. 


LATTICE CALCULATIONS AND RESULTS 


We study scattering using the same lattice action used to study the Hoyle State of uni, the structure of 

ll25l . and other nuclear Systems Il26l . The spatial lattice spacing is a = 1.97 fm and the temporal lattice spacing is at = 1.32 fm. 
Revisiting these calculations again in the future with different lattice spacings and going to higher order will provide a useful 
measure of systematic errors in lattice calculations of higher-body nuclear Systems EtII . We should mention that the dependence 
on lattice spacing has also been explored using model interactions of point-like alpha particles on the lattice Il28ll29]| . 

We perform projection Monte Carlo simulations with auxiliary fields to compute the matrices and on a 

periodic cubic lattice with length L = 16 fm. See Ref. OÜl for an overview of methods used in lattice effective field theory. We 
détermine from calculations with time steps and from calculations with L* + 1 time steps. The projection 

time T equals the product of Lt and the temporal lattice spacing Oj. For these calculations, a new algorithm is used to allow for 
Monte Carlo updates of the auxiliary fields as well as updates of the alpha cluster positions. 

We compute the radial adiabatic Hamiltonian using Eq. (|^ and extend the radial adiabatic Hamiltonian to a much larger 
volume with length 120 fm. This is done by computing matrix éléments of at large séparation from single-alpha lattice 

simulations and then including the Coulomb interaction between the otherwise non-interacting clusters. This process also allows 
us to define a “trivial” two-cluster Hamiltonian in which the two alpha clusters are non-interacting except for the infinite-range 
Coulomb interaction. 

With the radial adiabatic Hamiltonian defined in the large (120 fm)^ box, we extract the scattering phase shifts by imposing a 
hard spherical wall boundary at some chosen radius iîwaii and solving for the standing wave modes. In Fig.|^we show iS-wave 
radial functions for two different radial excitations (2S and 35”) at next-to-next-to-leading order in chiral effective field theory. 
We can extract the phase shift directly by fitting to the asymptotic behavior of the radial wave fonction as done in Ref. 113. 
However, we find it more efficient and accurate to extract the phase shifts from the energy différences between the spectrum of 
the adiabatic Hamiltonian and the trivial two-cluster Hamiltonian. This technique is discussed in detail in Ref. OTIl . 

Fig.j^shows the phase shifts for iS-wave scattering versus laboratory energy at next-to-leading order (NFO) in chiral effective 
field theory, next-to-next-to-leading order (NNFO), and the comparison with experimental data ll2Tlj24ll . The results at leading 
order (LO) do not include Coulomb effects and so hâve significantly different behavior near threshold and are not shown. We 
note that the NFO and NNFO phase shifts are quite similar, and both are in good agreement with experiment. The error bands 
are estimated from the spread in phase shifts versus the number of time steps in the projection Monte Carlo calculation, Lt- 

In Fig. I^we plot the Æ'-wave phase shifts at NNFO for Lt = A to Lt = 10. The error bars are computed using a jaekknife 
analysis of the stochastic errors of the Monte Carlo data. We find that the calculated ®Be ground State is bound at NNFO, though 
only a small fraction of an MeV away from threshold. For further comparison, we show in the inset next-to-leading-order results 
using halo effective field theory with point-like alpha particles lf3^ . We note that the norm matrix in Eq. 0 requires an even 
number of time steps in order to construct the inner product of dressed cluster States and so is not strictly defined for odd Lt- 
This explains some systematic différences between the even Lt and odd Lt results. However, they do appear to agréé for large 
Lt as expected. The data is not sufficiently précisé to allow for a multi-parameter extrapolation to the limit Lt ^ oo. So instead 
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FIG. 2. S-wave scattering wave functions at next-to-next-to-leading order versus radial distance. 



£L.b (MeV) 

FIG. 3. S-wave phase shifts at NLO and NNLO and comparison with experimental data 1214241 . 


the error band indicates the spread in results for different values of Lf. 

In Fig. I^we show phase shifts for U-wave scattering versus laboratory energy at NLO, NNLO, and the comparison with 
experimental data. We find that the 2+ résonance energy is pushed upwards by the Coulomb interactions and other higher-order 
corrections at NLO and NNLO. The results at NNLO are in fairly good agreement with the experimental results. 

In Fig. 1^ we plot the D-wave phase shifts at NNLO for Lt = 4 to Lt = 10. Again the error bars are computed using a 
jackknife analysis of the stochastic errors of the Monte Carlo data. In this case we see a stronger systematic différence between 
even Lt and odd Lt results, with the even Lt data approaching from above and the odd Lt data approaching from below. 


SUMMARY AND OUTLOOK 

In this letter we describe the first ab initia calculation of ^He + "‘He scattering. We use lattice effective field theory and the 
adiabatic projection method to compute phase shifts for ^-wave and Z?-wave scattering up to NNLO and find good agreement 
with experimental data. In order to perform these calculations, several new innovations were developed. This includes the use 
of spherical wave projections of the lattice initial States and a new algorithm that performs updates of both the auxiliary field 
configurations and alpha cluster positions. We plan to revisit these ^He + ^He calculations again in the future with different 
lattice spacings and going one order higher to NNNLO. These phase shifts provide useful benchmarks to assess systematic errors 
in calculations of higher-body nuclear Systems. 

Perhaps the most significant outcome of this study is that there is now a new method for scattering and reactions that has 
very favorable scaling with particle number. The computational effort needed for the Ai-body + 2 l 2 -body problem is roughly 
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5--wave at NNLO 



ÊLab (MeV) 

FIG. 4. 5-wave phase shifts at NNLO for different numbers of projection time steps Lt and comparison with experimental data 12111241 . In 
the inset we show next-to-leading-order calculations using halo effective field theory fÿ2\ . 



FIG. 5. O-wave phase shifts at NLO and NNLO and comparison with experimental data 12111241 . 


(^1 + ^ 2 )^ for light and medium-mass nuclei, and there is no restriction on the projectile being very light. Since the problem 
of sign oscillations is greatly suppressed for alpha-like nuclei lfÏ4l[^lMl . this approach appears to be a viable method to study 
important alpha processes involving heavier nuclei. 
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